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Abstract 

Dynamics and separation of mixed chiral microswimmers are numerically investigated in a chan¬ 
nel with regular arrays of rigid half-circle obstacles. For zero shear flow, transport behaviors are 
the same for different chiral particles: the average velocity decreases with increase of the rotational 
diffusion coefficient, the direction of the transport can be reversed by tuning the angular velocity, 
and there exists an optimal value of the packing fraction at which the average velocity takes its 
maximal value. However, when the shear flow is considered, different chiral particles show differ¬ 
ent behaviors. By suitably tailoring parameters, particles with different chiralities can move in 
different directions and be separated. In addition, we also proposed a space separation method 
by introducing a constant load, where counterclockwise and clockwise particles stay in different 
regions of the channel. 
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I. INTRODUCTION 


Active matter in biological and physical systems has been studied theoretically and 
experimentally|l-3|. Active particles or agents are assumed to have an internal propul¬ 
sion mechanism, which may use energy from an external source and transform it under 
non-equilibrium conditions into directed motion. Understanding this kind of motion can 
provide insight into out-of-equilibrium phenomena associated with biological examples such 
as bacteria, as well as with artificial microswimmers. These microswimmers may appear 
in several heldsj^, 3], e.g., to perform tasks in lab-on-a-chip devices, to localize pollutants 
in soils, and to deliver drugs within tissues. Compared with passive Brownian particles, 
active particles moving in confined structures could exhibit peculiar behaviors j-4 - 35], result¬ 
ing for example in spontaneous rectified transport^ -8], phase separation of self-propelled 
particles jib 111. tr app ing of particles in microwedge 12], depletion of elongated particles from 


low-shear regions 


131 ] . spira 


vortex formation in circular confinement 
16] . 


3 , 


and collective mo¬ 


tion in complex systems | 

Separation strategies of active particles based on their swimming properties 


20, are of 


utmost importance for various branches of science and engineering. Costanzo and coworkers 
[uj proposed a separation method for run-and-tumble particles in terms of their motility in a 
microchannel. The centrifugation can be also used to separate swimming cells having differ¬ 
ent motility 19] ■ Yang et. al. 20] proposed to use self-driven artificial microswimmers for the 
separation of binary mixtures of colloids. McCandlish and coworkers 21] studied the sponta¬ 
neous segregation of active and passive particles, freely swimming in a two-dimensional box. 
Asymmetric obstacles can also induce the separation of self-propelled particles with differ¬ 
ent motilities 22]. These separation strategies were mainly based on different motilities. 


However, in some situations, chirality separation is very important, for example often only 
one specific chirality is needed by the pharmaceutical and chemical industry. Therefore, a 
deeper understanding of chiral sorting mechanisms is of great fundamental and technologi¬ 
cal importance. Recently, Mijalkov and Volpe 23] demonstrated that chiral microswimmers 
can be sorted on the basis of their swimming properties by employing some simple static 
patterns in their environment. In asymmetrically patterned arrays, particles with the same 


radius but different chirality move in different directions 


24], 


In this paper, we numerically study the dynamics and separation of mixed chiral particles 
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in a channel with regular arrays of rigid half-circle obstacles. Two chiral separation methods 
are proposed: (A)Movement direction separation: in the presence of the shear flow, the 
counterclockwise particles move to the right while the clockwise particles move to the left. 
(B)Space separation: in the presence of a constant load, the counterclockwise particles stay 
in the upper region, while the clockwise particles stay in the bottom region. 


II. MODEL AND METHODS 


In the presence of an additional torque, the microswimmer tends to execute circular 
obits called as a chiral microswimmer. For example, when E. coli bacteria and spermatozoa 
undergo helicoidal motion near boundaries, which becomes two-dimensional chiral active 


Brownian motion 


25|. The kinetic of chiral active particles in conhned structures could 


exhibit peculiar behavior. In a channel, chiral mciroswimmers drift autonomously along 
a narrow channel under more general asymmetry conditions 26] and the diffusion of chi¬ 
ral mciroswimmers depends on both the strength and the chirality of the self-propulsion 
mechanism 27]. 



FIG. 1: (Color online) (a) Scheme of the separating device: two kinds of chiral self-propelled par¬ 
ticles moving in Couette flow with shear rate 7 . The half circle obstacles are regularly arrayed in 
the channel. Periodic boundary conditions are imposed in the x direction, and hard wall bound¬ 
aries in the y direction. The red and blue balls denote counterclockwise and clockwise particles, 
respectively. (b)The effective entropic potential V e ff( r x): the reduction of the coordinates (i. e. 
elimination of y coordinate by assuming equilibrium in the orthogonal direction) can involve the 
appearance of the entropic barrier. 


In this paper, we consider mixed chiral particles moving in Couette flow 


34] with walls 
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(the height H) in the ^-direction and periodic boundary conditions (the period L) in the 
x-direction shown in Fig. 1(a). Strictly, we do not consider microswimmers but active 
particles, because the hydrodynamic coupling between particle motion was ignored. The 
rigid half-circle obstacles with the diameter D are regularly arrayed in the channel. The 
active particles are modeled as soft disks with the diameter a and it is viewed as characterized 
by a unit vector n* = (cos 9 *, sin 9f) in the x — y plane, defining the direction of the self- 
propelled velocity. The effect of the shear flow enters the equations via the shear rate 7 . 
The forces acting particle i from the other particles and from the obstacles are respectively 
defined as F, = Ffe x + Ffe y = JT and G* = Gfe x + G^e y = JT G ty The dynamics of 
particle i is described by the following overdamped Langevin equations 


7 Ui + ^0 COS 9i + + G* + flood), 

(1) 

^ = v 0 sin 9i + y(Ff + Gf), 

(2) 

+ \Z^F>ofi{t), 

( 3 ) 


where /j is the mobility and r>o is the magnitude of the self-propelled velocity. Do is the 
rotational diffusion coefficient, flood, is a constant load along x direction. The Gaussian white 
noise £j(£) satisfies (£i(t)) = 0, and (£i{t)£j(s)) = 5ijS(t — s). is the angular velocity and the 
sign of determines the chirality of the particle. We define particles as the counterclockwise 
particles for positive and the clockwise particles for negative 

The particle-particle interaction and particle-obstacle interaction are assume to be of the 
linear spring form with the stiffness constant k 0 and k, respectively. F^ = k 0 (a — Vij)e r , 
if rij < a (F,j = 0 otherwise), and ry,- is the distance between particle i and j. G^ = 
k[\(a + D) — Tij\e r , if r t j < |(a + D) (Gy = 0 otherwise), and r t j is the distance between 
particle i and obstacle j ( for a contact with the flat side of a half circle, D — 0). 

Upon introducing characteristic length scale and time scale x — y — t — fik 0 t, 
Eqs. (flll2l[3l) can be rewritten in the dimensionless forms 


dxi 

dt 


IVi + ^0 COS 9i + Ff + Gi + flood-. 


dyj 

dt 


vo sin 9i + Ff + G \, 


( 4 ) 

( 5 ) 
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( 6 ) 




and flood = From now on, we will use only the dimensionless variables and shall omit 


the hat for all quantities occurring in the above equations. Note that the model with zero 


shear flow and zero constant load was also studied in the previous work 35], where self- 


propelled particles can present a vortex-type motion around convex symmetric obstacles 
even in the absence of hydrodynamics effects. 

The ratio between the area occupied by particles and the total available area can be 
defined as the packing fraction <f = Nna 2 /[A(LH — aS' q )] , where S 0 is the area covered by 
the obstacles and N is the total number of particles. The system is simulated in a box with 
the above described boundary conditions. In the asymptotic long-time regime, the average 
velocity of particles along x direction can be obtained from the following formula 



( 7 ) 


i=1 


For convenience, we define the scaled average velocity V s = v/vq through the paper. 

In order to facilitate discussion, we can obtain the effective entropic potential after 
the elimination of y coordinate by assuming equilibrium in the orthogonal direction and 
V e ff(x) (X —k,B ln^f, here k B is the Boltzmann constant and S(x) is transverse cross sec¬ 
tion and defined as 



( 8 ) 


and the effective potential is shown in Fig. 1(b). 

III. RESULTS AND DISCUSSION 

In this work, we propose two methods for chirality separation: (A) Movement direction 
separation by shear flow and (B)Space separation by a constant load. Unless otherwise 
noted, our simulations are under the parameter sets: D = 30, L = 40, H = 40, and k = 100. 
For the convenience of discussion, we focus on equimolar mixtures, where the number of the 
counterclockwise particles is equal to the number of the clockwise particles. 
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A. Chirality separation induced by shear flow 


In this section, we will systematically study the scaled average velocity V s in terms of Dg, 
|h 2 |, (j), and 7 in the presence of shear flow at fi oa d = 0 . 0 . 




FIG. 2: (Color online) Average velocity V s as a function of the rotation diffusion Dg at cf) = 0.1. 
(a) 7 = 0 and |f2| = 0. (b )7 = 0.004 and |0| = 0.02. 


Figure 2(a) shows the average velocity V s of nonchiral particles as a function of the 
rotation diffusion Dg without the shear flow (7 = 0.0). It is found that directed current 
occurs and the curves are almost the same for both clockwise (described by the blue line) 
and counterclockwise (described by the red line) particles. When Dg is very small, the self- 
propelled angle 6 almost does not change, V s tends to a saturate value. As Dg increases, 
V s decreases monotonically. For very large Dg, particles change their directions very fast. 
The self-propelled velocity acts as a zero mean white noise and the average velocity then 
tends to zero. However, the transport behavior becomes different when the shear flow is 
considered (shown in Fig. 2(b)). On increasing Dg from zero to 0.1, the average velocity 
of the clockwise particles increases monotonically from a negative value, while the average 
velocity of the counterclockwise particles decreases from a positive value. When Dg > 0.1, 
the average velocities are the same for two kinds of particles. Obviously, particles in shear 
flow can be separated by the chirality only for small values of Dg(e. g. Dg < 0.002). 
This behavior can be explained as follows. Due to the existence of the half circle, the 
counterclockwise and clockwise particles move easily in the up and down regimes of channel, 
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respectively. When the shear flow is considered, the counterclockwise (clockwise) particles 
feel the positive (negative) shear force on average. Therefore, the counterclockwise particles 
move to the right while the clockwise particles move to the left. 




FIG. 3: (Color online)Average velocity V s as a function of the angular velocity |f2| at (f> = 0.1 and 
Dg = 0.0005 . (a )7 = 0. (b)y = 0.004. 


The average velocity V s as a function of the angular velocity |fl| is shown in Fig. 3. For 
zero shear flow (7 = 0.0)(shown in Fig. 3(a)), the term vq cos0(t)(oc vq cos(|0|t)) in Eq. ([I]) 
can be seen as the external driving force and the effective potential is shown in Fig. 1(b). In 
the adiabatic limit |fi| —>■ 0 , the external force can be expressed by two opposite static forces 
Vo and — vq, yielding the mean velocity V s = ^[V(v 0 ) + V(—'i>o)]- In this case, it is easier for 
particles moving toward the slanted side than toward the steeper side (see Fig. 1(b)), so the 
average velocity is positive. On increasing |0|, due to the higher frequency, particles do not 
get enough time to cross the slanted entropic barrier (the right side) which is at a larger 
distance from the minima. Since the distance from the minima to the basin of attraction 
of next minima from the steeper side (the left side) is less than from the slanted side (the 
right side), hence in one period particles get enough time to climb the entropic barrier from 
the steeper side than from the slanted side, which resulting in a negative current. When 
| 0 | —>• 00 , the self-propelled angle changes very fast, particles will experience a zero averaged 
constant force, so V s tends to zero. When the shear flow is involved, the average velocity 
of counterclockwise particles is always positive and it increases with the angular velocity. 
However, for the clockwise particles, the average velocity is positive for small and large 
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values of |0|, negative for medium values of |fi|. Therefore, there exists an interval of |fi| 
where the counterclockwise particles move to the right while the clockwise particles move 
to the left. 




FIG. 4: (Color online) Average velocity V s as a function of the packing fraction 0 at Dg = 0.0005 
. (a )7 = 0.0 and |O| = 0.0. (b )7 = 0.004 and |fl| = 0.02. 


Figure 4(a) shows the average velocity V s as a function of the packing fraction 0 for zero 
shear flow. The interaction between particles can cause two results: (A)activating motion 
in analogy with thermal noise activated motion for a single stochastically driven ratchet, 
which facilitates the ratchet transport and (B) reducing the self-propelled driving, which 
blocks the ratchet transport. When the packing fraction increases from zero, the factor A 
first dominates the transport, so the average velocity increases with the packing fraction. 
However, when 0 > 0.6, the factor B becomes more important, the average velocity decreases 
with increasing 0. Therefore, there exists an optimal value of 0 at which the average velocity 
takes its maximal value. For the case of 7 = 0.004 and |O| = 0.02 (shown in Fig. 4(b)), 
we can find that when 0 < 0.4 the counterclockwise particles move to the right while the 
clockwise particles move to the left (from Fig. 3(b)). For this case, the average velocity of 
the counterclockwise particles is negative (positive) and its magnitude decreases (increases) 
with increasing 0 for 0 < 0.4(0 > 0.4). On increasing 0, the average velocity of the clockwise 
particles first increases, and then decreases. 

Figure 5(a) depicts the average velocity V s as a function of the shear rate 7 with half circle 
obstacles. It is found that the average velocity of the counterclockwise particles increases 








Y 


Y 


FIG. 5: (Color online)Average velocity V s as a function of the shear rate 7 . (a)With half circle 
obstacles (D = 30.0). (b)Without half circle obstacles (D = 0.0). The other parameters are 
\n\ = 0 . 02 , cj) = 0 . 1 , and D e = 0.0005. 


monotonously with the shear rate. However, on increasing 7 from zero to 0.006, the average 
velocity of the clockwise particles is negative and its magnitude first increases to its maxima, 
and then decreases to zero. When the shear rate is large enough(e. g. 7 > 0.006), the 
average velocities for both counterclockwise and clockwise particles are always positive. 
This is because the large shear rate dominates the transport and the chirality of the particle 
can be neglected. Therefore, chirality separation can be realized only for small shear rate. 

Note that both the shear flow and the obstacles are required for chirality separation in 
this setup. The shear flow alone does not suffice to drive the chiral species in opposite 
directions. For example, in the absence of the obstacles, the average velocity goes to zero 
for both counterclockwise and clockwise particles shown in Fig. 5(b). Obviously, without 
the obstacles, no net current occur even in the presence of the shear flow. 

Finally, we discuss the effect of the particle size on the transport when fixing the particle 
number(iV = 160) shown in Fig. 6 . When a —> 0, the average velocity tends to a constant, 
which shows that the ratchet effect still appear even for point particles. When a >(L~ 
.D)/ 2 . 0 , particles are too large to pass through the obstacles, thus the average velocity tends 
to zero. Therefore, there exists an optimal value of a at which the average velocity takes its 
maximal value. 
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FIG. 6 : (Color online) Average velocity V s as a function of the diameter a of particles at N = 160. 
The other parameters are |fl| = 0.02, 7 = 0.0, and Dq = 0.0005. 


B. Chirality separation driven induced by a constant load 

In the above section, we proposed a movement direction separation method, where the 
counterclockwise particles move to the right while the clockwise particles move to the left 
under the special parameters. Now we will propose a space separation method where the 
counterclockwise and clockwise particles stay in different regions of the channel. In this 
method, we add a constant force fi oa d along x direction and no the shear flow is needed 

(7 = 0). 

Figure 7(a) shows the configuration of particles at t — 0. The corresponding particle 
density in y direction is shown in Fig. 7(c). At t = 0, both counterclockwise and clockwise 
particles are randomly distributed in the channel. These two kinds of particles are completely 
mixed together. After t = 1000, the distribution of particles changes (shown in Figs. 7(b) 
and 7(d)): the counterclockwise particles are distributed in the upper region, while the 
clockwise particles stay in the bottom region. Therefore, we can separation these two kinds of 
particles by space distribution. In order to give the explanation of the separation mechanism, 
we plot the schematic diagram of the space separation in Fig. 8 . As shown in Fig. 8 , due 
to the existence of the load fi oa d, all particles move to the left on average. When particles 
meet the flat side of a half circle, different kinds of particles have different behaviors. The 
counterclockwise particles can easily pass across the half circle barrier through the upper 
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FIG. 7: (Color online)(a)Initial system configuration with uniform particle density at t = 0. (b)The 
same system at t = 1000. (c)Particle density p(y) in y direction at t = 0. (d)Particle density p(y) 
in y direction at t = 1000. The red (blue) circle denotes the counterclockwise (clockwise) particle. 
The other parameters are |fl| = 0.05, 4> = 0-2, Dg = 0.0005, and fi oa d = —0.1. 

region, while the clockwise particles pass across the barrier through the bottom region. 
Therefore, after a long time, the counterclockwise (clockwise) particles will be gathered in 
the upper (bottom) region. 

Note that the direction of the constant load is important for separating particles. The 
particle separation occurs only when the particles are driven towards the flat edges of the 
obstacles (e. g. / < 0 in the present setup). When particles are driven towards the rounded 
edges of the obstacles (e. g. / > 0), two kinds of particles have almost the same space 
distribution (shown in Fig. 9), thus mixed particles can not be separated. 
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FIG. 8: Schematic diagram of the space separation. The red (blue) circle denotes the counter¬ 
clockwise (clockwise) particle. 



y 

FIG. 9: Particle density p{y ) in y direction at t = 1000. The other parameters are |f2| = 0.05, 
cf) = 0.2, D e = 0.0005, and fi oad = 0.1. 

C. Chirality separation in 3D systems 

The separation mechanism for the case of two dimensional chiral particles performing 
circular motion can be extended to the case of three-dimensional (3D) chiral particles. The 
schematic diagram of the 3D setup is shown in Fig. 10, where mixed chiral particles moving 
in a channel with walls in the y and ^-directions and periodic boundary condition in the re¬ 
direction. The half-circle obstacle is replaced by the half-cylinder obstacle with the diameter 
D and the height L z in 3D setup. The active particles are modeled as soft spheres with the 
diameter a and the packing fraction is rewritten as (j) = Nna 3 /[6L Z (LH — 0.1257tD 2 )]. In 
the x — y plane, we define particles as the counterclockwise particles for positive D and the 
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clockwise particles for negative f2. To transfer the two-dimensional equations to the three 
dimensional case, the orientation vector n = (sin ipi cos 6 1 *, sin ipi sin 6i, cos ipi ) is used and the 
third component z is added to all quantities 36]. The other parameters are the same as that 
in the two-dimensional case. 



FIG. 10: (Color online)Schematic diagram of the 3D separating device: two kinds of chiral self- 


propelled particles moving in the channel with walls in the y and z-directions and periodic boundary 


condition in the x-direction. The half-cylinder obstacle is regularly arrayed in the channel. 




y 

FIG. 11: (Color online)Particle separation in the 3D setup. (a)Movement direction separation by 
the shear flow (corresponding to Fig. 5(a)). (b)Space separation by a constant load (corresponding 
to Fig. 7(d)). The red (blue) circle denotes the counterclockwise (clockwise) particle. The height 
of the half-cylinder L z = 4.0. The other parameters are the same as the corresponding 2D case. 
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The results for 3D separating devices are shown in Figs. 11(a) and 11(b). It is found 
that two separation methods are still valid when L z « L (e. g. Lz = 4.0). When L z —y 0, 
the 3D setup will reduce to 2D setup. From the extensive numerical simulations (not shown 
in the paper), the ratchet velocity V s decreases with increase of L z . From the experimental 
point of view, the relatively large ratchet velocity V s is necessary for particle separation. 
Therefore, the condition of L z « L must be considered when we design the 3D separating 
devices. 

From the above results, we can find that the optimal parameters for particle separation, 
D g < 0.002 and 0.01 < |D| < 0.1 in the dimensionless forms. From the experimental point 
of view, it would be interesting to estimate the real values of parameters. The characteristic 


17] are a = 2 pm, v 0 =10/im/s, 


parameters of Janus particles in a low Reynolds number regime 
and /iko =100 s _1 . From the scaling relations, one can obtain the real values of the rotational 
diffusion coefficient (D e < 0.2 rad 2 /s) and the angular velocity (1 rad/s< |fl| <10 rad/s). 
Obviously, it is possible to experimentally realize our separating mechanisms with these 
characteristic parameters. 

Note that our results are obtained from the projection of the half-cylinder in x — y plane. 
However, in 3D situation, a clockwise particle with its z vector pointing upward is identical 
to a counterclockwise particle with its z vector pointing down. Therefore, in order to improve 
the distinguishability of chiral swimming in 3D setup, a constant force F z is applied to the 
z direction. In the presence of F z , all particles accumulate at the upward (down) plane with 
z = L z (z — 0) and can be easily distinguished. 


IV. CONCLUDING REMARKS 

To conclude, we numerically studied the transport of mixed chiral particles moving in 
regular arrays of rigid half-circle obstacles. The effect of the shear ffow enters the equations 
via the shear rate 7 . For zero shear ffow, both counterclockwise and clockwise particles 
have the same transport behavior. The average velocity decreases with increase of the 
rotational diffusion coefficient. The direction of the transport can be reversed by tuning the 
angular velocity. There exists an optimal value of the packing fraction at which the average 
velocity takes its maximal value. However, when the shear flow is considered, transport 
behaviors for different chiral particles become different. By suitably tailoring the system 
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parameters (the rotational diffusion Dq, the angular velocity |0|, the packing fraction (j), 
and the shear rate 7 ), the counterclockwise particles move to the right while the clockwise 
particles move to the left. In addition, in order to obtain the space separation of chiral 
particles, we considered the constant load along the ^-direction, instead of the shear flow. 
In this case, the counterclockwise particles are distributed in the upper region, while the 
clockwise particles stay in the bottom region, thus, we can separate chiral particles. The 
models are extended to the case of 3D chiral particles and the corresponding separating 
methods are still valid when L z « L. We expect that our separation methods can provide 
insight into out-of-equilibrium phenomena and have the technological applications for active 
microswimmers. 
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